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Abstract 

Background: Nearly 40 years ago, Freeland and Janzen predicted that liver biotransformation enzymes dictated 
diet selection by herbivores. Despite decades of research on model species and humans, little is known about the 
biotransformation mechanisms used by mammalian herbivores to metabolize plant secondary compounds (PSCs). 
We investigated the independent evolution of PSC biotransformation mechanisms by capitalizing on a dramatic 
diet change event — the dietary inclusion of creosote bush {Larrea tridentata) — that occurred in the recent 
evolutionary history of two species of woodrats {Neotoma lepida and N. bryanti). 

Results: By comparing gene expression profiles of two populations of woodrats with evolutionary experience to 
creosote and one population naive to creosote, we identified genes either induced by a diet containing creosote PSCs 
or constitutively higher in populations with evolutionary experience of creosote. Although only one detoxification 
gene (an aldo-keto reductase) was induced by both experienced populations, these populations converged upon 
functionally equivalent strategies to biotransform the PSCs of creosote bush by constitutively expressing aldehyde and 
alcohol dehydrogenases. Cytochromes P450s, methyltransferases, glutathione S-transferases and sulfotransferases. The 
response of the naive woodrat population to creosote bush was indicative of extreme physiological stress. 

Conclusions: The hepatic detoxification system of mammals is notoriously complex, with hundreds of known 
biotransformation enzymes. The comparison herein of woodrat taxa that differ in evolutionary and ecological 
experience with toxins in creosote bush reveals convergence in the overall strategies used by independent species 
after a historical shift in diet. In addition, remarkably few genes seemed to be important in this dietary shift. The 
research lays the requisite groundwork for future studies of specific biotransformation pathways used by woodrats to 
metabolize the toxins in creosote and the evolution of diet switching in woodrats. On a larger level, this work advances 
our understanding of the mechanisms used by mammalian herbivores to process toxic diets and illustrates the 
importance of the selective relationship of PSCs in shaping herbivore diversity. 
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Background 

For herbivorous animals, food ingestion can have deleteri- 
ous consequences [1,2]. At every meal, herbivores confront 
potential toxins in the form of plant secondary compounds 
(PSCs). Some herbivores, e.g., the black swallowtail {Papilio 
polyxenes) and Stephens woodrat {Neotoma stephensi), 
have evolved the ability to specialize on plant species with 
high concentrations of PSCs, whereas other herbivores 
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must consume undefended plants or small quantities of 
many species of plants with differing PSCs to keep toxin 
intake below threshold doses [3-6]. 

In general, diet selection and diet breadth of herbivores 
is thought to be governed by the capacity of biotransform- 
ation (detoxification) enzymes in the liver [7,8]. Major 
evolutionary dietary shifts seem to be enabled by an 
underlying change in biotransformation enzymes [9,10]. 
Surprisingly few alterations among myriad biotransform- 
ation enzymes are required for a change in the capacity to 
metabolize a new toxin. Insect resistance to pesticides and 
mammalian resistance to rodenticides can occur through 
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the increased activity of a single biotransformation enzyme 
[9,11]. Currently, far more is known about the biotrans- 
formation mechanisms employed by insect herbivores in a 
diet shift than their mammalian counterparts [1,9,12]. 

The vast diversity of biotransformation enzymes pre- 
sents a challenge in determining which enzymes enable 
the processing of a particular diet or toxin. For example, 
more than 300 hepatic biotransformation enzymes have 
been described for the laboratory rat, Rattus norvegicus 
(www.reactome.org, [13]). Enzyme activity assays exist 
for far fewer than the known number of biotransform- 
ation enzymes and are not always specific to a gene 
product [14]. To overcome these challenges, we took a 
genomic approach and used microarray technology cus- 
tomized for our study species, to identify, on a more glo- 
bal scale, biotransformation gene expression patterns 
involved in the processing of PSCs. Gene expression is an 
important phenotypic character, subject to selection from 
the environment [14-16]. Many recent studies have con- 
cluded that differences in gene expression among groups 
are the result of disparate selective forces [17-19]. 

We investigated the independent evolution of bio- 
transformation mechanisms related to diet switching by 
capitalizing on a dramatic diet change event that oc- 
curred in the recent evolutionary history two species of 
woodrats {Neotoma lepida, N, bryanti). Both species 
underwent radical dietary changes due to floral shifts 
during the climate change event of the late Pleistocene 
and early Holocene (18,700-10,000 years ago). This event 
resulted in both of these woodrats independently switch- 
ing from diets of juniper and/or cactus, to one contain- 
ing creosote bush as it naturally invaded from Mexico 
and expanded to become the dominant shrub in the 
major deserts of the U.S. southwest [20,21]. Fossil re- 
cords of middens of both species indicate that N, bryanti 
came into contact with creosote 7,000 years before N, 
lepida [22,23]. 

The change in diet to creosote from previous plant 
species represents a marked change in PSC composition. 
Juniper contains high concentrations of numerous ter- 
penes as well as less abundant tannins [24,25]. The pri- 
mary PSC in cactus (oxalate) is only degraded by gut 
microbes, not liver enzymes [26]. In contrast, creosote 
leaves produce a complex resin containing >300 com- 
pounds, primarily aromatic ones, which account for 10- 
25% of the leaf dry weight [27]. The principal component 
of resin is nordihydroguaiaretic acid (NDGA), a phenolic 
lignan with detrimental effects when administered to la- 
boratory rats in doses regularly consumed by woodrats 
[28-30]. The distinct differences in the PSCs of juniper 
and cactus versus creosote suggest they are metabolized 
by different biotransformation pathways [14,31]. Com- 
parative data on enzyme activities and gene expression of 
populations of N, lepida support this contention. Neotoma 



lepida that feed on juniper have different enzyme activity 
and gene expression profiles compared to those that feed 
on creosote even when they are fed the same diet [32-34] . 
Indeed, populations of N, lepida in the Mojave desert have 
adapted to a diet of creosote as evidenced by their ability 
to ingest greater quantities of creosote compared to those 
from the Great Basin desert that feed on juniper and have 
no prior exposure to creosote [35]. Less is known about 
the ancestral diet of N, bryanti and no data exist with re- 
spect to its ability to metabolize PSCs. However, a recent 
comparative study on gene expression in N, lepida and N, 
bryanti on a non-toxic diet revealed similarities in expres- 
sion profiles in populations of both species that feed on 
creosote, relative to N, bryanti populations naive to creo- 
sote [36] . These results are consistent with potential con- 
vergence in biotransformation strategies of N, lepida and 
N, bryanti with respect to metabolism of a creosote diet. 

To advance our knowledge of the mechanisms involved 
in diet switching in mammals as well as biotransformation 
of PSCs in general, we addressed three questions. We 
asked which genes are induced (upregulated) by woodrats 
in response to the addition of creosote PSCs to the diet. 
Gene induction is a conventional pharmacological ap- 
proach used to identify potentially important enzymes in 
biotransformation of foreign compounds. Second, we asked 
whether experience mattered with respect to the biotrans- 
formation of creosote. We contrasted the responses of 
woodrats that historically and currently consume creosote 
to a population that is ecologically and evolutionarily naive 
to creosote compounds. Lastly, we asked whether there 
was convergence in the pathways used for biotransform- 
ation of creosote resin between N, bryanti and N, lepida 
that regularly fed on creosote bush by comparing the in- 
duced transcripts and also by contrasting their gene ex- 
pression profiles. 

Results 

Feeding trial 

All three populations maintained weight on both the 
control and creosote diet (Additional file 1: Table SI). 
However, there was a significant difference across popu- 
lations with the experienced N, bryanti being heavier 
than the other two populations. 

Microarray quality control 

Four of the 24 arrays failed to pass all 9 of Agilent's quality 
metrics. Three of these passed 8 of 9 metrics, and 1 passed 
7 of 9. All were kept in the analysis. On average, less than 
0.3% of the features across all arrays were flagged as non- 
uniform. The clustergram grouping individual woodrats 
by overall expression profile (all 6286 genes) grouped all 
the N, lepida in one monophyletic clade and each N, 
bryanti population in its own monophyletic clade, but did 
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not group individuals by diet within species designations 
(Additional file 2: Figure SI). 

Transcripts induced by creosote 

Both species and all populations altered gene expression 
in response to the creosote diet compared to the control 
The experienced N, bryanti had fewer induced (n = 7) 
and repressed transcripts (n = 7) than either the experi- 
enced N, lepida (induced n = 26, repressed n = 20) or the 
naive N, bryanti (induced n = 18, repressed n = 9); how- 
ever, the difference among the groups was not significant 
(X^ = 1.23 p = 0.54). There were no induced transcripts 
shared by all groups in this analysis. The experienced 
and naive populations of N, bryanti shared two induced 
transcripts, one of which was for a biotransformation 
gene (aldo-keto reductase 7 A3, Table 1). None of the 
treatment groups on creosote induced more than five 
different biotransformation genes (Table 1). The largest 
fold induction of a biotransformation transcript was ob- 
served in the naive N, bryanti, which induced sulfotrans- 
ferase 3A1 by 22x. All other biotransformation transcripts 
across all groups exhibited inductions of 7-fold or less. 

Differential expression between naVve and experienced 
wood rats 

There were significant differences in expression related 
to experience with creosote between N. bryanti popula- 
tions. Although the total number of transcripts with 
greater expression in both the experienced and naive 
was similar (199 in experienced vs. 183 naive, [Table 2 A 
with list of genes in Additional file 3: Table S2]), the 
experienced population expressed absolutely twice as 
many transcripts for biotransformation enzymes (Table 2A, 
= 14.72, df = 1, p < 0.001). Proportionally, 37% of all 
transcripts expressed to a greater degree by experienced 
N. bryanti were biotransformation related compared to 
only 18% of those in the naive N. bryanti. Many of the 
differentially regulated biotransformation transcripts 
coded for the same Genbank accession. The effective 
gene numbers, calculated from Shannons H, reflected 
the same pattern as the total biotransformation tran- 
script counts. For naive N. bryanti the effective gene 
number was 19.7 compared to an effective gene num- 
ber for experienced N. bryanti of 34.8. 

The degree to which upregulated genes were expressed 
was greater, overall, in the experienced N. bryanti (Table 3, 
Additional file 3: Table S2). The biotransformation tran- 
script with the greatest expression in naive N. bryanti 
and the only one with more than 10-fold higher expres- 
sion in naive compared to experienced animals, was a 
cytochrome P450 (2C65). In contrast, there were three bio- 
transformation transcripts in the experienced N. bryanti 
that exceeded 10-fold higher expression compared to 
naive N. bryanti. These transcripts are related to the 



functionalization of aldehydes (aldo-keto reductase 1C12) 
and glucuronic acid and glutathione conjugation path- 
ways (UDP glucuronosyltransferase 2B34, glutathione 
S -transferase, mu 7). 

Significant Gene Ontology (GO) terms (within biological 
process) were different in the naive versus experienced N. 
bryanti. Before comparison, terms across ontology tiers 
were filtered to include only those with significant, positive 
z-scores (>2) and ordered by the number of genes in each 
ontology term. GO terms are reported only if they include 
10 or more genes from the array. Naive N, bryanti had 
more than twice the number of GO terms that passed 
these criteria. The top terms related to responses to stim- 
uli; many of the other terms related to stress, damage and 
cell death (Table 4A). The experienced N, bryanti had 
fewer significant terms; the function of these related 
largely to metabolic processes (Table 4B). 

Differential expression between experienced woodrats of 
two species 

The experienced N, bryanti had increased expression of 
more transcripts than N, lepida (134 vs 109) as well as a 
greater degree to which those genes were expressed 
(Additional file 4: Table S3). This result was also true 
with respect to the biotransformation transcripts alone; 
N, bryanti had greater expression of 3x more transcripts 
compared to N, lepida on creosote and, in general, these 
genes were expressed to a much higher degree (Table 2B 
and Table 5). The effective gene numbers also reflect this 
pattern. The effective gene number with greater expres- 
sion in N, bryanti was 20.7, whereas the effective gene 
number for N, lepida was 14.4. 

In experienced N. bryanti, transcripts with the greatest ex- 
pression were related to conjugation of metabolites with 
glucuronic acid (different glucuronosyltransferases) followed 
by those related to conjugation with glutathione (gluta- 
thione S -transferases; Table 5). In N, lepida, the tran- 
scripts with the highest expression were those related to 
functionalization of aldehydes by aldehyde dehydrogenase 
(1 LI) and oxidation by Cytochrome P450 (2C65). Com- 
pared to N, bryanti, N, lepida also had greater expression 
of some of the same transcripts that were significantly in- 
duced by creosote feeding, i.e., aldehyde dehydrogenase 
(1 LI) glutathione S-transferase mu-1 (Tables 1 and 5). 

Significant GO terms overrepresented by experienced 
N, bryanti fed creosote were characterized by metabolic 
and catabolic processes as well as responses to stimuli 
(Table 6A). GO terms overrepresented by experienced 
N, lepida are characterized by responses to stimuli and 
signaling functions (Table 6B). 

Discussion 

Despite the vast knowledge of drug-metabolizing enzymes 
in humans and model species, the biotransformation 
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Table 1 Transcripts induced by creosote diet in three treatments: A. naive N. bryanti, B. experienced N. bryanti and 
C. experienced N. lepida 



Gene ID 


Gene description 


Ratio 


p-value 


Induced in N. bryanti ■ 


- naive 






NM_020565 


Sulfotransferase family 3A, member 1 


22.34 


0.024497169 


NM_021391 


Protein phosphatase 1, regulatory (inhibitor) subunit lA 


12.42 


0.005341544 


NM_1 73295 


UDP glucuronosyltransferase 2 family, polypeptide B17 


3.3 


0.01004607 


NM_027153 


Pirin, mRNA 


3.18 


0.049005275 


NM_1 73295 


UDP glucuronosyltransferase 2 family, polypeptide B17 


2.95 


0.005109673 


NM_007631 


Cyclin Dl 


2.75 


0.047341015 


NM_010145 


Epoxide hydrolase ^, microsomal 


2.51 


0.004864803 


NM_008761 


*FXYD domain-containing ion transport regulator 5 transcript variant 2 


2.51 


0.005930014 


NM_019144 


Acid phosphatase 5, tartrate resistant 


2.38 


0.030204487 


NM_0 16740 


SI 00 calcium binding protein Al 1 


2.37 


0.002381142 


NM_0 19693 


HLA-B-associated transcript lA 


2.31 


0.002135102 


NM_013215 


**Aldo-keto reductase family 7, member A3 


2.2 


0.048796617 


NM_009673 


Annexin A5 


2.11 


0.022303224 


NM_0 10664 


Keratin 18 


2.08 


0.034267738 


NM_028070 


AlkB, alkylation repair homolog 4 (E. coli) 


2.07 


0.026611663 


NM_0 13899 


Translocase of inner mitochondrial membrane 10 homolog (yeast) 


2.04 


0.001531745 


NM_0 13058 


Inhibitor of DNA binding 3 


2.03 


0.014931095 


NM_001111030 


Activin A receptor, type IC 


2.02 


0.015088922 


Induced in N. bryanti ■ 


- experienced 






NM_053346 


Neuritin 1 


3.99 


0.016579939 


NM_00 1164627 


Rho GTPase activating protein 8 transcript variant 1, MutualBestHitTo 


3.1 


0.0205116 


NM_013215 


**Aldo-keto reductase family 7, member A3 


2.93 


0.005358317 


NM_001109171 


Leucine rich repeat containing 20 


2.65 


0.001865456 


NM_008761 


*FXYD domain-containing ion transport regulator 5 transcript variant 2 


2.23 


0.023632655 


NM_027582 


RIKEN cDNA 4921521 F21 gene 


2.11 


0.02627791 1 


NM_001 106470 


Similar to KIAA1627 protein 


2.05 


0.027174423 


Induced in N. lepida - 


experienced 






NM_027406 


Aldehyde dehydrogenase 1 family, member LI 


6.9 


0.032217115 


NM_029662 


Major facilitator superfamily domain containing 2 


5.01 


0.009297156 


NM_027406 


Aldehyde dehydrogenase 1 family, member LI 


4.5 


0.023937762 


NM_001 014058 


Ubiquitin specific peptidase 18 


3.56 


0.02525724 


NM_001 184980 


Sulfotransferase family 2A DHEA-preferring member 5, SimilarTo 


3.43 


0.016254811 


NM_031004 


Smooth muscle alpha-actin 


3.37 


0.025794253 


NM_001 126273 


AlkB, alkylation repair homolog 2 (E. coli) 


3.26 


0.047075352 


NM_1 77200 


SV2 related protein homolog (rat)-like 


3.25 


0.007639427 


NM_010145 


Epoxide hydrolase 1, microsomal 


2.99 


0.013289403 


NM_1 78686 


Centrosomal protein 120 


2.92 


0.012963163 


NM_0 10358 


Glutathione S-transferase, mu 1 


2.8 


0.022522125 


NM_00 1168541 


Tsukushin transcript variant 1, MutualBestHitTo 


2.75 


0.007592024 


NM_022331 


Ubiquitin-like domain member 1 


2.64 


0.004966041 


NM_001 184980 


Sulfotransferase family 2A DHEA-preferring member 5, SimilarTo 


2.45 


0.036858859 


NM_031768 


Integrin, alpha E, epithelial-associated 


2.4 


0.035330416 
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Table 1 Transcripts induced by creosote diet in three treatments: A. naVve N. bryanti, B. experienced N. bryanti and 
C. experienced N. lepida (Continued) 



NM_ 


.198780 


Phosphoenolpyruvate carboxykinase 1 (soluble) 


2.4 


0.047233237 


NM_ 


_011393 


Solute carrier family 1, member 2, transcript variant 3 


2.37 


0.000661494 


NM_ 


J 45 123 


Cartilage acidic protein 1 


2.27 


0.028029577 


NM_ 


.133626 


Ribosome binding protein 1 


2.27 


0.042370631 


NM_ 


_029494 


RAB30, member RAS oncogene family 


2.19 


0.030712931 


NM_ 


.153392 


Tetratricopeptide repeat domain 39A, transcript variant 2 


2.15 


0.044515732 


NM_ 


.138953 


Elongation factor RNA polymerase II 2 


2.13 


0.017500695 


NM_ 


_022602 


Pim-3 oncogene 


2.06 


0.018759684 


NM_ 


.053433 


Flavin containing monooxygenase 3 


2.03 


0.01753571 


NM_ 


_028116 


Pygopus 1 


2.02 


0.002651983 


NM_ 


.021390 


Zinc finger protein Salll 


2.0 


0.006073039 



Bolded entries have known detoxification function. Asterisks indicate transcripts that are induced in both categories A and B. There are no shared induced 
transcripts between B and C. 



mechanisms used by mammalian herbivores to metabolize 
PSCs are largely unknown. Until recently, many studies 
were limited to analysis of one or a few biotransformation 
enzymes [32,37-40]. This study took advantage of micro- 
array technology customized for a unique study system to 
address this deficit in our understanding. We investigated, 
on a transcriptomic scale, the genes induced by a particu- 
lar suite of PSCs (creosote resin) and evaluated whether a 
common set of genes were expressed by herbivores with 
varying levels of evolutionary experience with these PSCs. 
We explored whether experienced herbivores have inde- 
pendently converged on regulation of a similar set of bio- 
transformation genes. We found that biotransformation 
enzyme expression does indeed vary with ecological and 
evolutionary experience with creosote and that independ- 
ent woodrat lineages employ, in part, similar strategies for 
successfully dealing with these shared PSCs. In addition, 
the results narrow the field from hundreds of possible bio- 
transformation genes to less than ten candidates. The 
work provides a testable framework for the changes in ex- 
pression of biotransformation enzymes that may have oc- 
curred as woodrats shifted from one toxic diet to another. 



Induced biotransformation genes 

We identified a narrow set of candidate genes relevant 
to the biotransformation of creosote. Surprisingly few 
biotransformation transcripts were induced by any of 
the groups fed creosote resin compared to the control 
diet. Of the hundreds of biotransformation enzymes, only 
four unique biotransformation transcripts were induced in 
the naive N. bryanti, six in the experienced N. lepida, and 
in the experienced N. bryanti, only a single biotransform- 
ation enzyme, an aldo-keto reductase (AKR7A3) was 
induced (Table 1). All of these transcripts encode for 
enzymes that act on substrates similar to the com- 
pounds present in creosote, particularly aromatic com- 
pounds [27,31]. These enzymes could function in tandem 
as a pathway to produce the final metabolite excreted in 
urine and/or feces [41]. For example, epoxide hydrolase 
acts on aromatic compounds such as naphthalene epoxide 
whose parent compound, napthalene, is present in creo- 
sote [41]. AKR enzymes, in turn, act on the metabolites 
produced by epoxide hydrolase, whereas glutathione S- 
transferases and UDP glucuronosyltransferases add conju- 
gates to the metabolites of AKR [42]. Thus, although few 



Table 2 Number of transcripts with significantly different expression in a comparison of (A) naVve and experienced 
N, bryanti on the creosote diet and (B) experienced N. bryanti and N, lepida on the creosote diet 



A. 



Treatment population 

Higher expression, N. bryanti, naive 
Higlier expression, N. bryanti, experienced 



Transcripts 



All 

183 
199 



Detox 
35 
76 



Chi-square 
11.92, df= 1, p< 0.001 



Treatment species 



Higher expression, N. bryanti, experienced 
Higher expression, N. lepida, experienced 



All 

134 
109 



Detox 
69 
23 



Chi-square 
= 14.72, df=l, p = 0.0001 



Chi square analyses compare the number of detoxification genes with higher expression to the overall number of genes with higher expression across the two 
experimental groups. 
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Table 3 Transcripts associated with detoxification function with significantly greater expression in (A) naive N, bryanti 
relative to experienced N. bryanti and (B.) experienced N. bryanti relative to naVve N. bryanti 



Gene ID Gene description 

Greater expression in naive N. bryanti 

Phase I - catalyze oxidation, reduction and hydrolysis reactions 
Alcohol dehydrogenases - oxidize alcohols 

NM_009626 C57BL/6 J alcohol dehydrogenase class 4 

Carboxylesterases - hydrolyze carboxylic acid esters 



Carboxylesterase 2 
Carboxylesterase 2 
Carboxylesterase 1 
Carboxylesterase 1 
Carboxylesterase 5, 



NM_1 45603 
NM_1 45603 
NM_021456 
NM_021456 
NM_1 72759 

Cytochromes P450 - oxidize wide range of organic substrates 

NM_028191 Cytochrome P450, family 2, subfamily c, polypeptide 65 

NM_007825 Cytochrome P450, family 7, subfamily b, polypeptide 1 

XM_219933 PREDICTED: P450 family 2 subfamily c polypeptide 79, SimilarTo 

NM_028191 Cytochrome P450, family 2, subfamily c, polypeptide 65 

NM_007825 Cytochrome P450, family 7, subfamily b, polypeptide 1, 

NM_019138 Cytochrome P450, family 7, subfamily b, polypeptide 1 

NM_019138 Cytochrome P450, family 7, subfamily b, polypeptide 1 

NM_010009 25-hydroxyvitamin D3 lalpha-hydroxylase 

Flavin containing monooxygenases - oxidize amines 
NM_018881 Flavin containing monooxygenase 2, mRNA 

NM_018881 Flavin containing monooxygenase 2, mRNA 

Phase II - catalyze transfer of conjugates to metabolites 
Acetyltransferases - transfer acetyl conjugate 

NM_001 161712 Glycine C-acetyltransferase transcript variant 2 SimilarTo 

UDP glucuronosyltransferases - transfer glucuronic acid conjugate 

NM_1 72881 UDP glucuronosyltransferase 2 family, polypeptide B35 

NM_1 53598 UDP glucuronosyltransferase 2 family, polypeptide B34 

NM_1 53598 UDP glucuronosyltransferase 2 family, polypeptide 834 

NM_001 029867 UDP glucuronosyltransferase 2 family, polypeptide B36 

Glutathione S-transferases - transfer glutathione conjugate 
NM_012577 Glutathione S-transferase pi 

Methyltransferases - transfer methyl group conjugate 
NM_0 16785 
NM_0 16785 
NM_0 16785 
NM_0 16785 
NM_022884 
NM_022884 
NM_026440 
NM_022884 
NM_1 77846 



Thiopurine methyltransferase 
Thiopurine methyltransferase 
Thiopurine methyltransferase 
Thiopurine methyltransferase 
Betaine-homocysteine methyltransferase 2 
Betaine-homocysteine methyltransferase 2 
RNA (guanine-7-) methyltransferase 
Betaine-homocysteine methyltransferase 2 
MKIAA0547 protein 



Ratio 



2.11 

2.52 
2.52 
2.45 
2.36 
2.12 

15.62 

4.48 

4.43 

3.79 

3.24 

2.46 

2.2 

2.12 

2.78 
2.66 



3.45 

3.93 
2.86 
2.6 
2.02 

2.67 

5.2 

5.04 

4.72 

4.3 

3.48 

2.15 

2.12 

2.11 

2.04 



p-value 



0.019591 

0.001284 
0.001213 
0.036796 
0.040035 
0.001463 

0.006602 
0.000458 
0.001591 
0.009243 
0.001260 
0.002067 
0.002338 
0.021096 

0.018212 
0.026835 



0.015595 

0.048799 
0.000686 
0.000899 
0.000218 

0.003827 

0.000140 
0.000190 
0.000223 
0.000073 
0.001257 
0.005312 
0.004311 
0.005705 
0.000481 
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Table 3 Transcripts associated with detoxification function with significantly greater expression in (A) naive N, bryanti 
relative to experienced N. bryanti and (B.) experienced N. bryanti relative to naVve N. bryanti (Continued) 



Sulfotransferases - tranfer sulfo group conjugate 

NM_001 1 84980 Sulfotransferase family 2A DHEA-preferring member 5, SimilarTo 

NM_001 1 01 586 Sulfotransferase family 2A DHEA-preferring member 5, SimilarTo 

NM_001 1 01 586 Sulfotransferase family 2A DHEA-preferring member 5, SimilarTo 

NM_001 1 01 534 CDNA clone IMAGE:905371 8 

Greater expression in experienced N. bryanti 

Phase I - catalyze oxidation, reduction and hydrolysis reactions 
Aldo-keto reductases - oxidize and reduce aldehydes and ketones 
NM_01 3777 Aldo-keto reductase family 1 , member CI 2 

NM_03061 1 Aldo-keto reductase family 1, member C6 

Aldehyde dehydrogenases - oxidize aldehydes 
NM_1 53543 Aldehyde dehydrogenase 1 family, member L2 

NM_031057 Aldehyde dehydrogenase 6 family, member Al 

Carboxylesterases - hydrolyze carboxylic acid esters 
NM_1 45603 Carboxylesterase 2 

NM_001 1 90346 Carboxylesterase 2 transcript variant 2, SimilarTo 

NM_1 45603 Carboxylesterase 2 

NM_001 044258 Similar to Carboxylesterase 2 

NM_001 044258 Similar to Carboxylesterase 2 

NM_001 1 90346 Carboxylesterase 2, transcript variant 2, SimilarTo 

NM_1 72759 Carboxylesterase 5 

Cytochromes P450 - oxidize wide range of organic substrates 



NM_023025 
NM_0 12730 
NM_012730 
NM_012730 
NM_153312 
NM_153312 
NM_022434 



Cytochrome P450, family 2, subfamily J, polypeptide 4 
Cytochrome P450, family 2, subfamily d, polypeptide 2 
Cytochrome P450, family 2, subfamily d, polypeptide 2 
Cytochrome P450, family 2, subfamily d, polypeptide 2 
Cytochrome P450, family 3, subfamily a, polypeptide 23/polypeptide 1 
Cytochrome P450, family 3, subfamily a, polypeptide 23/polypeptide 1 
Cytochrome P450, family 4, subfamily f, polypeptide 14 



Flavin containing monooxygenases - oxidize amines 



NM_001161765 
NM_001161765 
Miscellaneous Phase I 

NM_0 13626 
NM_00 1004086 



Flavin containing monooxygenase 5 transcript variant 1, SimilarTo 
Flavin containing monooxygenase 5 transcript variant 1, SimilarTo 



Peptidylglycine alpha-amidating monooxygenase 
Paraoxonase 3 

Phase II - catalyze transfer of conjugates to metabolites 
Acetyltransferases - transfer acetyl conjugate 

NM_053853 N-acetyltransferase 1 

NM_001 108278 Spermidine/spermine Nl-acetyltransferase family member 2 

NM_053853 N-acetyltransferase 1 

NM_053853 N-acetyltransferase 1 

NM_001 006995 Acetyl-Coenzyme A acetyltransferase 2 

NM_001 009657 Histone acetyltransferase 1 

NM_001 009657 Histone acetyltransferase 1 



4.75 
4.73 
4.69 
4.11 



22.98 
7.65 

3.61 
2.19 

6.91 
6.55 
6.21 
6.04 
5.68 
3.63 
2.9 

6.7 

3.38 

3.34 

3.25 

3.1 

2.92 

2.34 

2.46 
2.4 

2.23 
2.18 



7.05 
6.26 
5.39 
4.99 
3.73 
2.24 
2.04 



0.006821 
0.020208 
0.025323 
0.031920 



0.000139 
0.000074 

0.000215 
0.046613 

0.000264 
0.000143 
0.000575 
0.000027 
0.000028 
0.000532 
0.001251 

0.002710 
0.035947 
0.019677 
0.047435 
0.006788 
0.000334 
0.002319 

0.005855 
0.006221 

0.005322 
0.000331 



0.000017 
0.000021 
0.000007 
0.000010 
0.004937 
0.000395 
0.000004 



Malenke et at. BMC Ecology 2014, 14:23 
http://www.bionnedcentral.conn/1472-6785/14/23 



Page 8 of 16 



Table 3 Transcripts associated with detoxification function with significantly greater expression in (A) naive N, bryanti 
relative to experienced N. bryanti and (B.) experienced N. bryanti relative to naVve N. bryanti (Continued) 



UDP glucuronosyltransferases - transfer glucuronic acid conjugate 

NM_1 53598 UDP glucuronosyltransferase 2 family, polypeptide B34 

NM_1 53598 UDP glucuronosyltransferase 2 family, polypeptide B34 

NM_201642 UDP glucuronosyltransferase 1 family, polypeptide A6B 

NM_15281 1 UDP glucuronosyltransferase 2 family, polypeptide Bl 

NM_1 53598 UDP glucuronosyltransferase 2 family, polypeptide B34 

NM_001 191676 UDP glucuronosyltransferase 2 family polypeptide B34, SimilarTo 

NM_1 53598 UDP glucuronosyltransferase 2 family, polypeptide B34 

Glutathione S-transferases - transfer glutathione conjugate 

NM_026672 Glutathione S-transferase, mu 7 

NM_026672 Glutathione S-transferase, mu 7 

NM_1 33994 Glutathione S-transferase, theta 3 

NM_1 33994 Glutathione S-transferase, theta 3 

NM_008183 Glutathione S-transferase, mu 2 

NM_008183 Glutathione S-transferase, mu 2 

NM_001 024361 Similar to Glutathione S-transferase Al 

XM_001 47391 1 PREDICTED: Glutathione S-transferase Mu 2, SimilarTo 

NM_001 024304 Glutathione S-transferase mu 4 

NM_008183 Glutathione S-transferase, mu 2 

NM_001 009920 Glutathione S-transferase Yc2 subunit 

Methyltransferases - transfer methyl group conjugate 

NM_001 008299 RNA (guanine-7-) methyltransferase 

NM_022884 Betaine-homocysteine methyltransferase 2 

NM_001 1 06470 Similar to KIAAl 627 protein 

NM_010321 Glycine N-methyltransferase 

NM_1 73765 Aminoadipate-semialdehyde dehydrogenase 

NM_010321 Glycine N-methyltransferase 

NM_027334 Methyltransferase like 7A 



11.7 
7.22 
5.21 
3.37 
2.24 
2.18 
2.06 

10.47 

9.14 

8.12 

7.66 

4.77 

2.68 

2.54 

2.41 

2.11 

2.09 

2.04 

6.35 
3.45 
3.35 
3.26 
3.02 
2.84 
2.73 



0.000347 
0.013582 
0.000008 
0.043313 
0.004634 
0.005383 
0.007980 

0.038865 
0.034704 
0.003689 
0.002784 
0.000073 
0.015619 
0.040554 
0.000974 
0.026064 
0.001695 
0.031287 

0.000009 
0.001721 
0.002862 
0.000333 
0.000019 
0.000468 
0.002067 



All animals were fed a creosote diet. Within populations, results are organized by major detoxification enzyme class. Duplicates indicate the response of multiple 
probes for a given gene. 



transcripts were induced in woodrats fed creosote, those 
induced transcripts produce enzymes that act on sub- 
strates similar to those in resin. Furthermore, the candi- 
date genes identified could function in concert with one 
another in the biotransformation of PSCs in creosote. 

Few of the transcripts induced by the creosote diet 
were shared across woodrat groups. The conspecific 
populations shared induction of AKR7A3. The super- 
family of AKR enzymes act on a broad variety of sub- 
strates, particularly aldehydes and ketones, and the 7 A3 
isoform metabolizes aflatoxin Bl [42]. Epoxide hydrolase 
was induced by the naive N. bryanti and experienced N, 
lepiduy and also the experienced N, bryanti but in this 
case, at less than the 2-fold cut-off. This enzyme adds 
water to epoxides that could otherwise cause toxicity or 
mutation [31]. Epoxides are often formed during the 



biotransformation of aromatic hydrocarbons [31], which 
are common in creosote resin [27]. The shared induc- 
tion across all three groups of woodrats suggests epoxide 
hydrolase could be critical to the biotransformation of 
resin. Lastly, the naive N, bryanti and experienced N, 
lepida induced different sulfotransferases. This superfam- 
ily has high affinity to myriad substrates and metabolizes 
compounds (e.g., polyaromatic hydrocarbons) present in 
creosote resin [31]. Sulfotransferases biotransform xenobi- 
otics by the addition of a sulfur co-factor, generated from 
dietary cysteine. Because the availability of cysteine is 
often limited, sulfation capacity is often much lower than 
other conjugation pathways. It is possible that the woodrat 
diet contains ample cysteine or its precursors, serine and 
methionine. Consistent with this idea, both experienced 
populations upregulated methyltransferases involved in 
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Table 4 Gene ontology (GO) terms overrepresented In 
(A) naVve N, bryanti and (B) experienced N. bryanti on 
creosote diets 

Ontology List Gene set z-score 



"List" indicates tlie number of genes liiglily expressed witliin tlie term; "Gene set" 
indicates tlie total number of genes in that ontology included on the array. 



serine and methionine metabolism (Tables 3 and 5). Alter- 
natively, sulfation capacity may be greater in woodrats 
compared to other mammals. However, previous work 
suggests the capacity for sulfation is exceeded in both naive 
and experienced N. lepida at low levels of creosote inges- 
tion [43]. The role that sulfation plays in the biotransform- 
ation of creosote resin deserves further attention. 

Alterations in biotransformation mechanisms as a result 
of diet shifts 

The change by woodrats from an ancestral diet of juni- 
per and/or cactus to that of creosote likely required a 
different set of biotransformation enzymes to process 
the disparate suites of PSCs in these plants. The naive 
and experienced N, bryanti had increased expression of 
about the same number of genes on creosote compared 
to the control diet. However, within that pool of genes, 
the naive N, bryanti expressed fewer transcripts with 
biotransformation functions than the experienced ani- 
mals. A similar pattern has been observed elsewhere; 
woodrats feeding on a novel diet (i.e., animals naive to 
the diet) expressed relatively few biotransformation tran- 
scripts and more transcripts related to cellular function 
compared to feeding on their native diet [33]. The ani- 
mals on the novel diet were thought to be unable to 
marshal the appropriate biotransformation response, and 
the upregulation of genes with diverse cellular functions 
may have been to prevent or respond to the physiological 
consequences of the PSCs. Indeed, the functional analyses 
of the gene ontologies for differentially expressed genes in 
N, bryanti are consistent with this hypothesis (Table 4). 
The transcriptomes of naive woodrats were enriched in 
GO classes related to considerable physiological stresses 
(e.g., death, response to wounding). In contrast, those 
enriched in experienced N, bryanti were related to the 
processing of toxins. 

Functional convergence 

There is evidence for functional convergence between 
the two species with respect to the transcripts induced 
on a creosote diet. As discussed above, the transcript 
for epoxide hydrolase was the only biotransformation 
transcript significantly induced by both species; how- 
ever, it was expressed at levels below the 2-fold cut-off 
used in the initial analysis. The top biotransformation 
transcript induced by N, lepida (aldehyde dehydrogen- 
ase) and the only biotransformation transcript induced 
by N, bryanti (aldo-keto reductase) both encode for en- 
zymes that metabolize similar substrates, particularly 
aldehydes. Creosote resin contains at least 300 com- 
pounds, including aldehydes and compounds that may 
have aldehyde functional groups after biotransform- 
ation by other enzymes [27,44]. Aldehydes are reactive 
compounds that cause considerable cellular damage. 



A. naYve N. bryanti 



response to stimulus 


61 


1241 


3.66 


response to chemical stimulus 


30 


533 


3.13 


system development 


28 


587 


2.06 


response to stress 


27 


525 


2.44 


immune system process 


22 


265 


4.62 


apoptosis 


20 


350 


2.56 


cell death 


20 


367 


2.34 


death 


20 


370 


2.3 


programmed cell death 


20 


353 


2.52 


regulation of multicellular organismal process 


20 


328 


2.88 


regulation of apoptosis 


17 


287 


2.51 


regulation of cell death 


17 


297 


2.35 


regulation of programmed cell death 


17 


289 


2.47 


cellular response to chemical stimulus 


15 


222 


2.89 


immune response 


15 


140 


4.92 


regulation of developmental process 


15 


244 


2.5 


response to external stimulus 


13 


189 


2.75 


regulation of multicellular organismal 
development 


12 


206 


2.02 


defense response 


11 


147 


2.83 


positive regulation of developmental process 


11 


113 


3.82 


response to wounding 


11 


137 


3.09 


hemopoiesis 


10 


103 


3.63 


hemopoietic or lymphoid organ development 


10 


106 


3.52 


immune system development 


10 


115 


3.23 


negative regulation of apoptosis 


10 


145 


2.41 


negative regulation of cell death 


10 


153 


2.23 


negative regulation of programmed cell death 


10 


145 


2.41 


B. experienced N, bryanti 








metabolic process 


96 


2711 


2.39 


small molecule metabolic process 


35 


782 


2.51 


response to chemical stimulus 


25 


533 


2.3 


cellular ketone metabolic process 


23 


332 


4.23 


carboxylic acid metabolic process 


22 


321 


4.07 


organic acid metabolic process 


22 


328 


3.96 


oxoacid metabolic process 


22 


321 


4.07 


amine metabolic process 


13 


203 


2.81 


cellular amine metabolic process 


12 


181 


2.82 


cellular amino acid metabolic process 


11 


152 


3.02 


monocarboxylic acid metabolic process 


10 


157 


2.43 
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Table 5 Transcripts associated with detoxification function with significantly greater expression in (A.) experienced 
N, bryanti relative to N. lepida and (B.) experienced N. lepida relative to N. bryanti 



Gene ID Gene Description 

Greater expression in experienced N. bryanti 

Phase I - catalyze oxidation, reduction & hydrolysis reactions 
Aldehyde dehydrogenases - oxidize aldehydes 

NM_1 7871 3 Aldeliyde deliydrogenase 8 family, member Al 

Alcohol dehydrogenases - oxidize alcohols 

NM_017270 Alcoliol deliydrogenase 4 (class II), pi polypeptide 

NM_017270 Alcohol dehydrogenase 4 (class II), pi polypeptide 

NM_017270 Alcohol dehydrogenase 4 (class II), pi polypeptide 

NM_017270 Alcohol dehydrogenase 4 (class II), pi polypeptide 

Carboxylesterases - hydrolyze carboxylic acid esters 

NM_001 1 90346 Carboxylesterase 2 transcript variant 2, SimilarTo 

XR_033674 PREDICTED: similar to Carboxylesterase 2 

NM_001 190346 Carboxylesterase 2 transcript variant 2. SimilarTo 

NM_001 1 90346 Carboxylesterase 2 transcript variant 2. SimilarTo 

Cytochromes P450 - oxidize wide range of organic substrates 

NM_008898 P450 (cytochrome) oxidoreductase 

NM_1 47206 Cytochrome P450, family 3, subfamily a, polypeptide 9 

NM_1 47206 Cytochrome P450, family 3, subfamily a, polypeptide 9 

NM_153312 Cytochrome P450, family 3, subfamily a, polypeptide 23/polypeptide 1 

Flavin containing monooxygenases - oxidize amines 

NM_008030 Flavin containing monooxygenase 3 

Superoxide dismutases - dismutase superoxide to oxygen and water 

NM_017050 Superoxide dismutase 1 

Phase II - catalyze transfer of conjugates to metabolites 
UDP glucuronosyltransferases - transfer glucuronic acid conjugate 



Ratio 



2.22 

3.03 
2.82 
2.32 
2.26 

2.61 
2.28 
2.24 
2.04 

2.36 
2.18 
2.09 
2.08 

2.52 

2.17 



NM_ 


J 52811 


UDP g 


ucuronosyltransferase 2 family, polypeptide Bl 


26.84 


NM_ 


J 52811 


UDP g 


ucuronosyltransferase 2 family, polypeptide Bl 


26.2 


NM_ 


.153598 


UDP g 


ucuronosyltransferase 2 family, polypeptide B34 


13.03 


NM_ 


.153598 


UDP g 


ucuronosyltransferase 2 family, polypeptide B34 


10.82 


NM_ 


.153598 


UDP g 


ucuronosyltransferase 2 family, polypeptide B34 


8.52 


NM_ 


.152811 


UDP g 


ucuronosyltransferase 2 family, polypeptide Bl 


7.33 


NM_ 


.152811 


UDP g 


ucuronosyltransferase 2 family, polypeptide Bl 


7.05 


NM_ 


.153598 


UDP g 


ucuronosyltransferase 2 family, polypeptide B34 


4.58 


NM_ 


.173295 


UDP g 


ucuronosyltransferase 2 family, polypeptide B17 


4.23 


NM_ 


.173295 


UDP g 


ucuronosyltransferase 2 family, polypeptide B17 


4.16 


NM_ 


.173295 


UDP g 


ucuronosyltransferase 2 family, polypeptide B17 


3.82 


NM_ 


.153598 


UDP g 


ucuronosyltransferase 2 family, polypeptide B34 


3.79 


NM_ 


.153598 


UDP g 


ucuronosyltransferase 2 family, polypeptide B34 


3.2 


NM_ 


.153598 


UDP g 


ucuronosyltransferase 2 family, polypeptide B34 


3.13 


NM_ 


.001191676 


UDP g 


ucuronosyltransferase 2 polypeptide B34, SimilarTo 


2.8 


NM_ 


.153598 


UDP g 


ucuronosyltransferase 2 family, polypeptide B34 


2.64 


NM_ 


.001191676 


UDP g 


ucuronosyltransferase 2 family polypeptide B34 


2.52 


NM_ 


.201642 


UDP g 


ucuronosyltransferase 1 family, polypeptide A6B 


2.34 



p-value 



0.010400 

0.005872 
0.005802 
0.000857 
0.002058 

0.003223 
0.000942 
0.002730 
0.021685 

0.013012 
0.021999 
0.023839 
0.014854 

0.004188 

0.007301 



0.000202 
0.000230 
0.000721 
0.001924 
0.000934 
0.000471 
0.001293 
0.003637 
0.004670 
0.004114 
0.006675 
0.003846 
0.000542 
0.000852 
0.003254 
0.001950 
0.010507 
0.009381 
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Table 5 Transcripts associated with detoxification function with significantly greater expression in (A.) experienced 
N. bryanti relative to N. lepida and (B.) experienced N. lepida relative to N. bryanti (Continued) 



NM_009467 


UDP glucuronosyltransferase 2 family, polypeptide B5 


2.23 


0.003903 


NM_201642 


UDP glucuronosyltransferase 1 family, polypeptide A6B 


2.23 


0.005852 


Glutathione S-transferases - transfer glutathione conjugate 






NM_0 12796 


Glutathione S-transferase, theta 2 


5.9 


0.002555 


NM_0 12796 


Glutathione S-transferase, theta 2 


4.93 


0.003108 


NM_0 12796 


Glutathione S-transferase, theta 2 


4.79 


0.003690 


NM_00 1024361 


Similar to Glutathione S-transferase Al 


3.98 


0.005397 


NM_00 1024361 


Similar to Glutathione S-transferase Al 


3.58 


0.008732 


NM_00 1024361 


Similar to Glutathione S-transferase Al 


3.52 


0.005938 


NM_00 1024361 


Similar to Glutathione S-transferase Al 


3.25 


0.014456 


NM_0 12796 


Glutathione S-transferase, theta 2 


2.89 


0.032248 


NM_00 1024361 


Similar to Glutathione S-transferase Al 


2.66 


0.008524 


NM_00 1024361 


Similar to Glutathione S-transferase Al 


2.39 


0.010509 


NM_008183 


Glutathione S-transferase, mu 2 


2.34 


0.047800 


NM_00 1009920 


Glutathione S-transferase Yc2 subunit 


2.32 


0.018875 


NM_00 1009920 


Glutathione S-transferase Yc2 subunit 


2.1 


0.025697 


NM_001 077353 


Glutathione S-transferase, alpha 3, transcript variant 2 


2.08 


0.048535 


Methyltransferases 


- transfer methyl group conjugate 






NM_009349 


Indolethylamine N-methyltransferase 


10.85 


0.000877 


NM_009349 


Indolethylamine N-methyltransferase 


8.38 


0.002392 


NM_009349 


Indolethylamine N-methyltransferase 


6.52 


0.009050 


NM_009349 


Indolethylamine N-methyltransferase 


6.47 


0.007757 


NM_080462 


Histamine N-methyltransferase 


2.13 


0.004838 


XM_223974 


Methyltransferase 1 1 domain containing 1 


2.04 


0.005653 


NM_1 72687 


Coenzyme Q3 homolog, methyltransferase (yeast) 


2.03 


0.008959 


Sulfotransferases - 


tranfer sulfo group conjugate 






NM_020565 


Sulfotransferase family 3A, member 1 


4.17 


0.006575 


NM_020565 


Sulfotransferase family 3A, member 1 


3.81 


0.001338 


NM_020565 


Sulfotransferase family 3A, member 1 


3.52 


0.007965 


NM_0 18805 


Heparan sulfate (glucosamine) 3-0-sulfotransferase 3B1 


3.22 


0.001553 


NM_020565 


Sulfotransferase family 3A, member 1 


3.2 


0.000616 


NM_020565 


Sulfotransferase family 3A, member 1 


2.83 


0.004497 


NM_031834 


Sulfotransferase family, cytosolic, lA, phenol-preferring, member 1 


2.79 


0.001314 


NM_020565 


Sulfotransferase family 3A, member 1 


2.77 


0.006085 


NM_0 18805 


Heparan sulfate (glucosamine) 3-0-sulfotransferase 381 


2.75 


0.000909 


NM_020565 


Sulfotransferase family 3A, member 1 


2.53 


0.000740 


NM_031834 


Sulfotransferase family, cytosolic, lA, phenol-preferring, member 1 


2.26 


0.003965 


NM_031834 


Sulfotransferase family, cytosolic, lA, phenol-preferring, member 1 


2.05 


0.004712 


NM_020565 


Sulfotransferase family 3A, member 1 


2.04 


0.024154 



B. Greater expression in experienced N. lepida 

Phase I - catalyze oxidation, reduction & hydrolysis reactions 
Aldehyde dehydrogenases - oxidize aldehydes 

NM_027406 Aldehyde dehydrogenase 1 family, member LI 9.62 0.030847 

NM_027406 Aldehyde dehydrogenase 1 family, member LI 4.1 1 0.045209 
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Table 5 Transcripts associated with detoxification function with significantly greater expression in (A.) experienced 
N, bryanti relative to N. lepida and (B.) experienced N. lepida relative to N. bryanti (Continued) 



Alcohol dehydrogenases - oxidize alcohols 






NM_007410 


Alcohol dehydrogenase 5 (class III), chi polypeptide 


2.14 


0.025440 


Cytochromes P450 


- oxidize wide range of organic substrates 






NM_028191 


Cytochrome P450, family 2, subfamily c, polypeptide 65 


8.87 


0.017817 


NM_007820 


Cytochrome P450, family 3, subfamily a, polypeptide 16 


3.69 


0.006413 


NM_019138 


Cytochrome P450, family 7, subfamily b, polypeptide 1 


2.68 


0.015580 


NM_007815 


Cytochrome P450, family 2, subfamily c, polypeptide 29 


2.53 


0.028674 


NM_023025 


Cytochrome P450, family 2, subfamily J, polypeptide 4 


2.45 


0.004241 


NM_0 10009 


25-hydroxyvitamin D3 lalpha-hydroxylase 

Phase II - catalyze transfer of conjugates to metabolites 


2.37 


0.030007 


Glutathione S-transferases - transfer glutathione conjugate 






NM_010358 


Glutathione S-transferase, mu 1* 


2.99 


0.021783 


NM_010358 


Glutathione S-transferase, mu 1* 


2.94 


0.009552 


NM_017013 


Glutathione S-transferase A2 


2.38 


0.007870 


NM_017013 


Glutathione S-transferase A2 


2.34 


0.010162 


NM_017013 


Glutathione S-transferase A2 


2.29 


0.010172 


Methyltransferases 


- transfer methyl group conjugate 






NM_030241 


SET domain containing (lysine methyltransferase) 8 


2.19 


0.001281 


NM_025907 


Methyltransferase like 6 


2.05 


0.019380 


NM_0 16668 


Betaine-homocysteine methyltransferase 


2.04 


0.013408 


Sulfotransferases - 


tranfer sulfo group conjugate 






NM_001101534 


CDNA clone IMAGE:9053718 


2.21 


0.007280 


NM_027928 


Carbohydrate (chondroitin 4) sulfotransferase 13, MutualBestHitTo 


2.07 


0.006258 


NM_027928 


Carbohydrate (chondroitin 4) sulfotransferase 13, MutualBestHitTo 


2.04 


0.006510 


NM_027928 


Carbohydrate (chondroitin 4) sulfotransferase 13, MutualBestHitTo 


2.01 


0.009812 



Within species, results are organized by major detoxification enzyme classes. Duplicates Indicate the response of multiple probes for a given gene. 



Inadequacies in the biotransformation of aldehydes are 
connected to a number of human diseases [45,46]. 
Given the potential for damage, aldehydes in creosote 
are likely a significant selective force resulting in the 
convergence on the induction of transcripts related to 
aldehyde metabolism by both species. 

Further demonstration that the two experienced woo- 
drats exhibit a parallel response to the ingestion of creo- 
sote resin is witnessed in their overall gene expression 
patterns. Many of the biotransformation transcripts that 
were more highly expressed in the experienced N. 
bryanti have analogous functions with those that were 
more highly expressed in N, lepida (Table 5). For ex- 
ample, both species have higher expression of different 
isoforms of aldehyde dehydrogenases (ALDH8A1 versus 
ALDHILI). A similar pattern is present for alcohol dehy- 
drogenases. Cytochromes P450, methyltransferases, gluta- 
thione S-transferases, and sulfotransferases. One notable 
exception to this pattern is that N, lepida did not have a 
corollary in the UDP-glucuronosyltransferases (UGTs), a 
superfamily of biotransformation enzymes responsible for 



glucuronidation. Neotoma bryanti had higher expression 
of at least 6 different UGTs, one of which was expressed 
26 fold higher. The UGTs are considered some of the 
most versatile of the conjugation enzymes due to their 
wide range of substrates [31]. Previous studies reported 
that N, lepida fed creosote increased glucuronidation 
thereby demonstrating it is indeed an important pathway 
in N, lepida [32,43]. It is possible that both species use 
the same UGTs for creosote biotransformation but that N, 
bryanti utilizes this pathway to a greater extent resulting 
in the higher gene expression values observed in this 
study. 

The last line of evidence for convergence stems from 
GO results. The experienced populations fed creosote 
share three ontologies (response to stimulus, oxidation- 
reduction process, lipid metabolic process), which all re- 
late to biotransformation (Table 6). 

Constitutive differences 

While many biotransformation enzymes are induced 
when an animal is exposed to xenobiotics, some are 
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Table 6 Gene ontology (GO) terms overrepresented In 
(A) N, bryanti and (B) N. lepida on a creosote diet 

Ontology List Gene set z-score 



metabolic process 65 2711 2.58 

response to stimulus 36 1241 2.73 

small molecule metabolic process 30 782 4.09 

cellular response to stimulus 26 882 2.28 

catabolic process 24 626 3.57 

cellular catabolic process 21 550 3.28 

response to chemical stimulus 18 533 2.44 

oxidation-reduction process 17 375 3.68 

cellular ketone metabolic process 15 332 3.42 

carboxylic acid metabolic process 14 321 3.15 

organic acid metabolic process 14 328 3.06 

oxoacid metabolic process 14 321 3.15 

response to organic substance 14 348 2.82 

small molecule catabolic process 13 282 3.25 

lipid metabolic process 12 316 2.38 

cellular response to chemical stimulus 10 222 2.74 

response to endogenous stimulus 10 186 3.37 
B. experienced N. lepida 

response to stimulus 35 1241 2.13 

Signaling 25 764 2.52 

signal transduction 24 689 2.79 

regulation of response to stimulus 16 430 2.49 

intracellular signal transduction 15 347 3.03 

oxidation-reduction process 15 375 2.7 

regulation of signal transduction 15 327 3.28 

regulation of signaling 15 363 2.84 

lipid metabolic process 13 316 2.61 

"List" indicates the number of genes highly expressed within the term; "Gene set" 
indicates the total number of genes in that ontology included on the array. 

expressed constitutively at high levels [31]. It is difficult 
to identify such baseline differences in expression using 
microarrays. For example, a recent study on an herbivor- 
ous insect that specializes on glucosinolate-rich plants 
did not identify any induced biotransformation tran- 
scripts when the insects were fed glucosinolates com- 
pared to a control diet [47]. Constitutive expression 
could explain this outcome. In our study, direct compar- 
isons of gene expression between populations fed creo- 
sote provide insight into how animals differ in baseline 
biotransformation, especially when coupled with the re- 
sults from the induction study. The experienced N. 
bryanti induce very few transcripts, only one of which 
has known biotransformation function, relative to both 
the naive N, bryanti and the experienced N, lepida. 
However, direct comparisons of experienced N, bryanti 



with both the naive conspecific and experienced congen- 
eric when all are fed creosote, revealed many more dif- 
ferences in expression (Table 2 A and B). Overall, the 
experienced N, bryanti had greater expression of signifi- 
cantly more biotransformation enzymes. This pattern 
coupled with the paucity of induced transcripts suggests 
that the experienced N, bryanti may be constitutively ex- 
pressing biotransformation enzymes at a higher level. It 
is possible that N, bryantis longer historic exposure to 
creosote caused an increase in baseline production of 
biotransformation enzymes such that they continually 
express the appropriate combination of enzymes to effi- 
ciently detoxify creosote PSCs [22,23]. Moreover, the expe- 
rienced N. lepida GO terms include many for regulation 
and signal transduction, suggesting processes in flux, 
whereas these terms are absent from the experienced N. 
bryanti results (Table 6). 

Conclusions 

Herbivores and the plants on which they feed represent 
one of the most common interactions in nature, yet we 
know relatively little about the mechanisms that herbi- 
vores, particularly mammalian ones, employ to over- 
come plant toxins [1,48]. The application of recently 
developed genome-based approaches is enabling re- 
searchers to more thoroughly investigate the responses 
of herbivores to plant defensive compounds. This work 
capitalized on the availability of a recently developed 
microarray to identify the genes induced by mammalian 
herbivores when feeding on PSCs. The results lay the 
requisite groundwork for future functional and evolu- 
tionary studies of the genes involved in the metabolism 
of creosote toxins and the evolution of diet switching in 
woodrats. 

Methods 

Trapping and feeding trial 

We trapped two species, N, lepida (desert woodrat) and N, 
bryanti (Bryants woodrat). Neotoma lepida were trapped 
at Lytle Ranch Preserve (lat: 37.117514, long: -114.009661, 
Washington Co., UT, USA). Neotoma bryanti were trapped 
near Palm Springs (lat: 33.679616 long: -116.362018, 
Riverside Co., CA, USA). These two populations were 
chosen because both the desert woodrat and Bryants 
woodrat at these sites include creosote bush (L tridentata) 
as a primary component of their diet [35] (unpublished ob- 
servations by MDD and JRM). With respect to creosote in- 
gestion, we refer to these populations as "experienced," 
both in their individual histories as well as their evolution- 
ary histories. In preliminary trials, experienced animals 
from both species were capable of ingesting high levels of 
creosote resin (8% by dry mass) added to a powdered 
rabbit chow base (Teklad formula 2031). In addition, we 
trapped N, bryanti at Ronald W. Caspers Wilderness Park 



A. experienced N. bryanti 

metabolic process 65 2711 

response to stimulus 36 1241 

small molecule metabolic process 30 782 

cellular response to stimulus 26 882 

catabolic process 24 626 

cellular catabolic process 21 550 

response to chemical stimulus 18 533 

oxidation-reduction process 17 375 

cellular ketone metabolic process 15 332 

carboxylic acid metabolic process 14 321 

organic acid metabolic process 14 328 

oxoacid metabolic process 14 321 

response to organic substance 14 348 

small molecule catabolic process 13 282 

lipid metabolic process 12 316 

cellular response to chemical stimulus 10 222 

response to endogenous stimulus 10 186 

B. experienced N. lepida 

response to stimulus 35 1241 

Signaling 25 764 

signal transduction 24 689 

regulation of response to stimulus 16 430 

intracellular signal transduction 15 347 

oxidation-reduction process 15 375 

regulation of signal transduction 15 327 

regulation of signaling 15 363 

lipid metabolic process 13 316 
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(lat: 33.53367, long: -117.54965, Orange Co., CA, USA). 
This population has no evolutionary or ecological exposure 
to creosote bush, and is considered "naive" to creosote 
PSCs. In preliminary trials, these animals did not maintain 
body mass on diets of 4% creosote resin. All experimental 
procedures were approved by the University of Utahs Insti- 
tutional Animal Care and Use Committee (#07-02015). 

Eight woodrats from each of the three groups (experi- 
enced N, lepida, experienced N, bryantU naive N, bryanti) 
were divided into two treatments of four individuals (n = 2 
males, 2 females). One treatment was fed powdered rabbit 
chow amended with resin from creosote bush; the other 
group was fed a control diet of the same powdered rabbit 
chow with no additions. Resin was extracted and creosote 
diets prepared as in [34]. 

Animals given the resin diet were fed a gradually in- 
creasing amount of resin over an eight-day period to 
allow induction of biotransformation enzymes [31]. This 
group received powdered chow with 0% resin for three 
days, 1% creosote resin (dw/dw) for two days, followed 
by 2% creosote resin for three days. The final 2% con- 
centration was selected to accommodate the naive N, 
bryanti since they were unable to tolerate more resin 
without excessive mass loss (>10% initial weight). Ani- 
mals in the control treatments received 0% resin diet for 
all eight days. Animals had water ad libitum, and diet 
was presented daily. Animals were weighed daily and 
body masses were compared within and between groups 
using a repeated measures ANOVA with time as the re- 
peated variable and population and diet as factors. 

Animals were dispatched using CO2 asphyxiation; liver 
tissue was preserved in RNAlater (Sigma) and archived 
at -80°C. RNA was extracted (RNAqueous) for the 
microarray analysis. 

Microarray analysis 

Liver samples were analyzed using a custom microarray, 
built from the hepatic transcriptome of a single desert 
woodrat [36]. Two primary probe groups were designed 
from the transcriptome. The target contigs for the first 
group (n = 943) were woodrat contigs whose annotation 
matched a list of biotransformation genes extracted from 
Norway rat arrays that we had previously used in studies 
with woodrats [34,49]. These probes were all labeled 
"biotransformation." The target contigs in the second 
group consisted of contigs from the woodrat transcriptome 
that had >75% sequence identity with a rodent BLAST 
match in the region of highest homology (n = 5355). All 
probes were screened for quality and replicated on the 
array. 

The RNA quality and concentration of each sample was 
evaluated, and the Agilent One-Color Quick Amp Labeling 
Kit was used to generate fluorescently labeled cRNA for 
hybridizations. Additional details on microarray processing 



methods are available in [36]. Microarray features were 
extracted using Agilent Feature Extraction software ver- 
sion 10.5.1.1. All control spots, non-uniform spots and 
population outlier spots were removed from the data- 
set; intensity values were log^ transformed using Agi- 
lentPilter, software designed to simplify the processing 
of Agilent data. Data from duplicate probes were com- 
bined, but different probes with the same annotation 
were maintained separately since it is possible that the 
original contigs were from different isoforms despite 
the shared match to a rodent gene. Consequently, many 
of the resulting lists of differentially expressed genes have 
multiple seemingly identical entries. For most analyses, we 
treat these probes as targeting unique genes; in a few spe- 
cified analyses, we have combined or discarded duplicate 
probes to evaluate the results as conservatively as possible. 
The performance of the custom Neotoma expression array 
was previously validated [36]. 

Species and diet comparisons 

To compare gene expression across diet and species, the 
data for the Neotoma probes was batch uploaded to Gene- 
sifter 3.7. Prior to comparing expression profiles, the 
consistency of transcriptome response was evaluated by 
comparing overall gene expression profiles across individ- 
uals. After normalizing the intensity data, individuals were 
clustered by gene expression using all woodrat- derived 
probes (n = 6286). Clustering parameters were distance: 
correlation, linkage:average, row centered:by genes. 

We created a total of five pairwise comparisons in 
Genesifter. The specific nature of these comparisons 
are described in the subsequent paragraphs. In all com- 
parisons, the quality requirement was set to 1, and all 
individuals from both treatments were required to pass. 
T-tests were performed with alpha = 0.05, and resulting 
lists of significantly differentially expressed genes/probes 
were ordered by fold change with a fold change threshold 
of greater or equal to 2. This approach is less conservative 
than implementing a statistical control for the false dis- 
covery rate (e.g., Benjamini and Hochberg correction 
[50]), but fold change ranking is more consistent [51]. 

To identif)^ induction and repression of genes with re- 
spect to ingestion of resin, we compared expression on 
the creosote diet to the control diet for each of the three 
groups (experienced N, lepida, experienced N, bryanti, 
naive N. bryanti). The lists of significantly differently 
expressed genes were ordered by direction and fold 
change. Transcripts that were upregulated on the creosote 
diet were considered induced by creosote and may be in- 
dicative of enzymes important in the biotransformation of 
creosote. Transcripts that were downregulated were con- 
sidered repressed. To test for conserved or convergent re- 
sponses, the lists of induced genes were compared for 
entries present in more than one woodrat group. 



Malenke et at. BMC Ecology 2014, 14:23 
http://www.bionnedcentral.conn/1472-6785/14/23 



Page 15 of 16 



In the fourth comparison, naive and experienced N. 
bryanti fed the creosote diet, were compared to identify 
contrasts in expression with respect to prior experience. 
Using the gene lists generated from these comparisons, 
we compared the overall number of transcripts with 
greater expression in each group as well as the subset 
of these transcripts characterized by having biotrans- 
formation function. The standard gene lists treat probe 
entries with identical Genbank accession numbers as 
independent listings. To more conservatively compare the 
relative response of each group to the creosote diet, we 
quantified the number of unique Genbank accessions in 
the biotransformation gene lists and, then calculated the 
Shannons H Diversity Index from each list. To calculate 
Shannons H, each unique Genbank accession was treated 
as a "species" and the number of times it appeared in 
the significantly upregulated list was treated as its 
"abundance". The values for Shannons H were then used 
to calculate an effective number of upregulated genes i.e., 
the number of "species" expected in comparable commu- 
nities given their Shannons H. This index is used to 
contextualize diversity indices across multiple communi- 
ties [52]. The effective species value (or effective gene 
value, in this case) is calculated as exp (H). 

To examine function more broadly, we used Gene 
Ontology (GO). GO is a bioinformatics tool that groups 
individual genes by the function of their products. GO al- 
lows the organization of many individual gene results into 
fewer functionally-defined categories. There are 3 major 
classifications: cellular component, molecular function 
and biological process. Within each of these three, func- 
tion is categorized more and more specifically. Within our 
GO results, z-scores were used to determine significantly 
important associations. Ontology terms with z-scores 
of >2 are considered to be significantly overrepresented 
in the results; that is, the genes in that ontology are 
more likely to be upregulated than expected by chance, 
given the number features on the array within that par- 
ticular GO term. Within the biological process classifi- 
cation, two gene ontology lists were generated for 
transcripts with greater expression in naive N, bryanti 
and in experienced N, bryanti. Lastly, experienced ani- 
mals from both species fed the creosote diet were com- 
pared to illuminate contrasts in expression between 
species with similar experience. This contrast was con- 
ducted in the same fashion as the naive versus experi- 
enced N, bryanti. 

Additional files 



Additional file 1: Table SI. Animal masses before and after feeding trial. 

Additional file 2: Figure SI. Cluster analysis of individual woodrat 
samples labeled by treatment. Description: Branching pattern determined 
by overall gene expression patterns. Individual animals are labelled with 



species and experience with creosote bush in the wild; creosote diet in 
laboratory trial is indicated by bolded text. 

Additional file 3: Table S2. Significantly differently expressed transcripts 
from comparison of naive versus experienced N. bryanti. Description: (A.) 
Transcripts with greater expression by naive N. bryanti; (B.) transcripts with 
greater expression by experienced N. bryanti. All animals were fed a creosote 
diet. Duplicates indicate the response of multiple probes for a given gene. 

Additional file 4: Table S3. Significantly differently expressed 
transcripts from comparison of experience N. byranti versus experienced 
N. lepida. Description: (A.) Transcripts with greater expression by N. 
bryanti; (B.) transcripts with greater expression by N. lepida. All animals 
were fed a creosote diet. Duplicates indicate the response of multiple 
probes for a given gene. 
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